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Abstract. 

Low frequency spin fluctuation dynamics in paramagnetic spinel L1V2O4, a rare 
3d-electron heavy fermion system, is investigated. A parametrized self-consistent 
renormalization (SCR) theory of the dominant AFM spin fluctuations is developed and 
applied to describe temperature and pressure dependences of the low-T nuclear spin- 
lattice relaxation rate 1/Ti in this material. The experimental data for 1/Ti available 
down to ~ IK are well reproduced by the SCR theory, showing the development of 
AFM spin fluctuations as the paramagnetic metal approaches a magnetic instability 
under the applied pressure. The low-T upturn of 1/T\T detected below 0.6 K under 
the highest applied pressure of 4.74 GPa is explained as the nuclear spin relaxation 
effect due to the spin freezing of magnetic defects unavoidably present in the measured 
sample of L1V2O4. 
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1. Introduction 

The metallic spinel LiV204 has attracted much attention since a heavy fermion 
behaviour in this material was reported [1, 2, 3]. Despite continuous activity in last 
years, there is currently no consensus on the mechanism for formation of heavy fermion 
quasiparticles in LiV 2 04, and the issue is still under debate [4, 5]. 

At low temperatures, T < 30 K, the spin system of LiV 2 4 exhibits pronounced 
short-range antiferromagnetic (AFM) correlations [6, 7, 8], but no long-range magnetic 
ordering was detected at any measured temperatures. The geometrical frustration of 
the pyrochlore lattice of vanadium ions (in the mixed valence state V 3 ' 5+ ) is likely to 
be a crucial aspect of the problem. The frustration may suppress at any T a long- 
range ordering of strongly correlated itinerant electrons, but instead, the system is 
placed near to a magnetic instability. The emergence of largely degenerate low lying 
spin excitations in the ground state of LiV 2 4 is expected to be responsible for low- 
T properties of this material, including its heavy fermion behaviour. This appealing 
picture has been developed in detail in previous work [9, 10]. It can be examined 
by considering experimental results obtained by different techniques, like the nuclear 
magnetic resonance (NMR) and the inelastic neutron scattering (INS), probing low 
frequency spin fluctuations. 

The 7 Li-NMR studies of the spin fluctuation dynamics in LiV204 were reported 
in a series of papers [11, 12, 13, 14, 15]. For high temperatures, T > 60 K, the NMR 
relaxation data have been successfully explained [12] in terms of a V local moment 
formalism, while in the low-T region, T < 30 K, an approach based on the itinerant and 
strongly correlated electrons in the paramagnetic metal LiV 2 4 is more appropriate. 
Low-T measurements down to ~1K on samples of high purity under ambient pressure 
reveal a nearly constant value of the Knight shift K and linear T dependence of the 
nuclear spin-lattice relaxation rate T^ 1 as for normal metals, but with a very high 
value of (TxT)" 1 . An estimate [12] of the Korringa relation K 2 T{I '/ 'S = R, where 
S = hrfl/ ' AitkB^ni with 7 e and 7„ being the electronic and nuclear gyromagnetic ratios 
respectively, and R ~ 0.5 is less than unity, indicates [16, 17] that presumably AFM 
fluctuations are dominant in the spin-lattice relaxation at low T. 

The NMR measurements on high purity samples of LiV 2 04 under pressure up to ~ 5 
GPa were also reported by K. Fujiwara et al [18, 19]. For T < 10 K, the value of (TxT) -1 
which gives information on the q averaged dynamical spin susceptibility becomes larger 
on applying higher pressure and grows with decreasing temperature. At the same time, 
the Knight shift K probing only the static unuform susceptibility x (q = 0) was found 
to be nearly temperature independent and insensitive to the pressure above 2 GPa. 
These results were suggested [19] to be indicative of an increase under applying pressure 
of AFM spin correlations at some momenta q^O and their enhanced dominance over 
those at q = 0. 

If temperature is sufficiently low, the NMR properties of LiV 2 4 are strongly 
affected by a small amount of magnetic defects, n de f ect < 1 mol % [13, 14, 15]. The 
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relaxation of the longitudinal nuclear magnetization versus time is no more a single 
exponential one but described by a stretched exponential function with the characteristic 
relaxation rate l/T-j* showing a peak at some temperature T pea k ~ IK. Such a behaviour 
was proved [15] to originate from the spin freezing of magnetic defects below T pea k- 
With decreasing ndefect, the peak position of 1/T* is apparently shifted to lower T. 
Remarkably, a proper model analysis of the low-T NMR data obtained in powder 
samples of LiV 2 4 with varying ndefect < 1 m ol % has shown [15] that relaxation 
effects due to inhomogeneously distributed magnetic defects and homogeneous spin 
fluctuations inherent to magnetically pure LiV 2 4 are separable and thus can be 
examined independently. It is worth noting that in a single crystal of LiV204 containing 
a small amount of magnetic impurities or crystal defects, a somewhat different behaviour 
was observed [15], possibly because of a lack of the separability of the relaxation effects 
mentioned. 

A scenario explaining the considerable influence of a weak disorder on the 1ow-c<j 
spin dynamics detected by the low-T NMR measurements on LiV 2 4 was proposed 
by Johnston et al [14], following a more general consideration developed in [20]. As 
already noted, a critical aspect of the problem is the emergence in the ground state 
of pure LiV 2 4 of a large number of low lying spin excitations, implying a proximity 
of the system to a magnetic instability. Inhomogeneously distributed magnetic defects 
may locally lift the degeneracy of low lying spin excitations and cause their partial 
condensation, thus giving rise to a strong change of spin dynamics at sufficiently low 
temperatures. A large degeneracy of strongly enhanced and slow spin fluctuations 
in pure LiV 2 4 was confirmed by a combined analysis of the low temperture INS 
data [6, 7, 8] and the complementary calculations [9] of the dynamic spin susceptibility 
x(q, oj). In these calculations performed first at T = 0, the actual electronic band 
structure of LiV 2 4 obtained in the local-density approximation is used and effects 
of strong electron correlations are treated in the random phase approximation. As 
an extension for finite T, the self-consistent renormalization (SCR) theory [21] of 
spin fluctuations was proved [10] to be a helpful tool in explaining the temperature 
renormalization of the Iow-cj spin fluctuation dynamics in LiV 2 4 derived from INS 
measurements [6, 7, 8]. 

In the present study, the parametrized SCR theory is applied to elucidate the 
main features of the temperature and pressure dependences of the spin-lattice relaxation 
rate observed in the low-T NMR measurements on LiV 2 4 . As known [21], the SCR 
theory offers a phenomenological description for spin fluctuations in nearly ferro- or 
antiferromagnetic itinerant electron systems by taking into account effects of mode- 
mode coupling between spin fluctuations either at q = or q ^ 0, respectively; the 
latter case is applicable to the paramagnetic spinel LiV 2 4 . For details we refer to 
our recent work [10], where the basic equation of the SCR theory is solved numerically 
and the results are compared with INS data for LiV 2 4 . There, the values of empirical 
parameters entering the SCR theory are estimated to provide the best overall coincidence 
between the theory and INS experiment. 
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The outline of the paper is as follows. In section 2, an expression for the spin- 
lattice relaxation rate T^ 1 is derived in terms of the SCR theory and, first, an evolution 
of (TiT) -1 down to ~ 1 K with increasing pressure is examined. Next, since the 
temperature behaviour of (T{T)~ l down to much lower temperature ~ 60mK is available 
only at the highest applied pressure of 4.74 GPa, these particular data deserve a special 
attention. We argue, contrary to what K. Fujiwara et al suggested [19], that the 
upturn of (TiT) -1 detected below 0.6 K is not entirely due to homogeneous critical 
AFM spin fluctuations, but more likely is a signature of an additional nuclear spin 
relaxation mechanism, probably due to the spin freezing of magnetic defects. Summary 
and concluding remarks can be found in section 3. 

2. SCR theory for relaxation rate 1/Ti in LiV 2 4 
2.1. Background 

The nuclear spin-lattice relaxation rate due to electronic spin fluctuations is generally 
given by 

J_ 27gfc B T \- . 2 Imx(q, u n ) m 

where A q is a q dependent effective hyperfine coupling; in our calculations, the resonance 
frequency uo n will be taken in the limit uo n — > 0. In (1), the q summation is over the 
Brillouin zone (BZ) of the fee lattice inherent to the pyrochlore lattice of V atoms in 
the spinel structure L1V2O4; x(q, cj) is the dynamic spin susceptibility calculated per 
primitive cell (four V atoms) in units of (g/i^) 2 . Calculations suggest [9] a rather peculiar 
q-dependence of x (q, The resulting model for a distribution in q space of dominant 
spin fluctuations was checked [10] to provide a firm ground to describe the low-T INS 
measurements for LiV 2 4 . Main features of the model are discussed and used below to 
calculate 1/Ti. 

In the low-T limit, the paramagnetic state of LiV 2 04 is characterized by strongly 
enhanced and slow spin fluctuations occupying a large region in q space around the 
surface of a mean radius |q| ~ Q c ~ 0.6 A -1 , called the critical surface [9] (see Figure 
1). These strongly degenerate 1ow-cj AFM spin fluctuations dominate over those at 
smaller q. Such a rather peculiar distribution in q space of the dominant AFM spin 
fluctuations is a consequence of geometrical frustration of the pyrochlore lattice of V 
atoms, which can be traced back [9] to the underlying electronic band structure and 
the many sheet Fermi surface of the metallic spinel LiV 2 04. With increasing T (up 
to 60 K), the AFM fluctuations get suppressed while those at the BZ center remain 
nearly T independent. As noted in [10], the observed [6, 7, 8] warming shift of the 
low-cu integrated INS intensity (from |q| ~ Q c at T — > toward low q values at higher 
temperatures) does not require any significant temperature renormalization of x (q> ^) 
at small q. 
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2.2. 1/71 T obtained from SCR theory 

For the pure LiV204, two main contributions to the spin-lattice relaxation rate can be 
written as 



As discussed later, at T — > 0, the second contribution from AFM spin fluctuations at 
|q| ~ Qc is larger than the first one coming from the small q spin fluctuations. For 
finite but low T, because of a comparatively small variation with both temperature 
and pressure of the Knight shift, and hence, of the intrinsic uniform susceptibility 
X (q = 0, uj = 0), one expects for the low-u; and small q spin fluctuations much weaker 
temperature and pressure dependences than those observed [6, 7, 8] for AFM spin 
fluctuations. Therefore, in the subsequent analysis the contribution (1/TiT) is 
assumed to be a constant and treated below as an adjustable parameter of the fit 
procedure. In this approximation, we relate the experimentally observed temperature 
and pressure dependence of 1/T±T merely to that of dominant AFM fluctuations around 
the Q c surface. 

Let us choose a wave vector Q c with the end point lying on the Q c surface, and 
consider its vicinity, q' = Q c + q, (|q| <C |Q C |), as depicted in Figure 1. The inverse 
dynamic spin susceptibility can now be expanded as 

m^rm^ +AW2+Bi « ±)2 -' c "' (3) 

where g" and q -1 are components of q parallel and perpendicular to the normal Q C /|Q C | 
to the Q c surface. Since the dynamic spin susceptibility is considered around q = 0, 
only the leading term of C independent of q is taken in the expansion (3). This allows 
us to avoid an increase of the number of parameters in describing the dynamic spin 
susceptibility. 

In (3), A, B and C are empirical parameters that, together with \ (Qc5 T = 0), can 
be estimated [10] from INS and magnetic measurement data. The parameters A, B and 
C are usually taken to be T-independent in the low-T region where the SCR theory 
works well. Strong anisotropy in q space of AFM spin fluctuations, i.e. b = B/A <^ 1, 
is assumed and verified [10]. It is also reasonable using a spherical approximation to 
the Q c surface, implying that x(Qc',T) does not depend on the direction of Q c , i.e., 

x(Q c ;T) = x(Q c ;T). 

To be close to the standard notation of the SCR theory [21], we introduce, instead 
of A and C, the following parameters (here, h = ks = Q^b = 1 ): 

where q B is the effective radius of the BZ boundary given in terms of the lattice primitive 
cell volume v as qs = (67r 2 /vq) 1 ^ . Next, the reduced inverse susceptibility at |q| = Q c 
is defined as 

y 1 (T) = 2TMQ c ;Ty (5) 
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Figure 1. A cross-section in q space of the critical Q c surface (|Q C | — 0.6A -1 ) is shown 
schematically by the solid line. Strongly enhanced antiferromagnetic spin fluctuations 
dominating low-T properties of L1V2O4 are located in a close vicinity (of width 
Aq < 0.2A -1 ) to the Q c surface. The meaning of decomposition q' = Q c + q|| + q± 
for an arbitrary wave vector q' near the Q c surface is explained in the text. 



In this notation, one obtains for Imx (q, oj n ) / u) n in the limit u n — > and up to a 
constant factor the following expression 

Imx (Q c + q, co n ; T) 1 1 



Wn TnT A f. fm\ , /■_!!/_ \2 , , / , /_ x2l 2 



VQ (T) + (qW/q B Y + b(^/q B y 



(6) 



When performing the integration over q in (1), two dimensionless cutoff parameters, 
x c = (Qc I 1b) 2 and z c = (qc/qs), are introduced. For the former, the requirement x c < 1 
is sufficient. As the latter cutoff we take z c ~ 1/2, which distinguishes the region of the 
dominant AFM spin fluctuations from that near the BZ center. Near the Q c surface, 
the hyperfine coupling is assumed to be a constant Aq c . Then the resulting expression 
for (l/TiT)^ Qc reads as 

1 \ 3l 2 n h\A Q f (Q c \ 2 1 / 1 



T i T ) q ~ Qc *k B T T A \q B ) te c \v^(T) tan yftfiF) 

tan' 1 - ^ : 1, (7) 



y/y Q (T)+bx c ' y/y Q (T) + bx 
including again the proper dimensional constants. 

2. 3. Basic equation of SCR theory and empirical parameters 

The reduced inverse susceptibility yg (t), where t = T/T , obeys the following integral 
equation [10]: 







with 



ri(u) = lnT(M) - ( u- - j \nu + u- ^ln2vr, (9) 
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where V (u) is the gamma function. 

The present SCR theory includes a set of five parameters which are now denoted 
as dq (0), Ta, To, qq and bx c . The parameters T A and To characterize, at T — > 0, the 
momentum and frequency spread of the dominant AFM spin fluctuations, qq is the 
effective mode-mode coupling constant and bx c (<^ 1) describes a large anisotropy of the 
spin fluctuation dispersion (3) in q-space. From a fit to INS data, we obtained [10] the 
following estimates 

T ~ 60K, T A ~ 220K, g Q = 0.16, bx c = 0.01, (10) 

and y Q (0) ~ 0.044. 

The parameter yg (0) = [2TaX {Qc',T = 0)]" 1 is a measure of distance from the 
magnetic instability. Following earlier studies [22, 23], we consider yg (0) to be the only 
pressure dependent parameter, assuming that yg (0) — > as the system approaches 
a quantum critical point with increasing pressure. In case of ambient pressure, both 
the empirical parameter yg (0) ~ 0.044 and the solution of equation (8) are taken the 
same as in [10], thus providing the agreement with INS data. Here, two remaining 
fit parameters for 1/TT under ambient pressure are found to be (1/TiT) ~ 0.55 
(sec _1 K _1 ) and \Aq c \ ~ 5 kG. For NMR data obtained under the applied pressure, 
the only adjustable parameter yg (0) is estimated by solving first the equation (8) and 
inserting the solution yg (T) into (7). 

2.4- Comparison to experimental data 

The experimental data [19] for (1/TiT) measured for different applied pressure together 
with theoretical curves fitting these data are depicted in Figure 2. A good agreement 
between the theoretical results and the available experimental data is found at least down 
to 1 K. First we note that away from the instability, i.e., yg (0) > 0, the SCR theory 
predicts the Korringa behaviour (1/TiT = constant) in the low-T limit. On applying 
higher pressure (or, as yg (0) — > 0), 1/TiT becomes larger, while a temperature range 
near T =0, where the Korringa relation holds, shrinks and tends to zero at yg (0) — > 0. 
These results of the SCR theory are in accordance with the temperature and pressure 
dependence of (1/TiT) observed in L1V2O4. 

The SCR theory fails, however, in giving a quantitative description of the low-T 
upturn of 1/TT detected below 0.6 K under the highest applied pressure of 4.74 GPa. 
For instance, as seen from Figure 2, the parameter yg (0) =0.018 fits well in describing 
the experiment down to ~ IK, while for 0< yg (0) <0.018 a strong deviation in the whole 
temperature range is found (see, for instance, the upper theoretical curve corresponding 
to yg (0) =0.008). It is worth emphasizing that just at the quantum critical point, 
Dq (0) = 0> the equation (8) leads to the power law (1/TT) ~ j n a g reem ent with 

the earlier SCR theory predictions [22, 23] for the quantum critical behaviour of (1/TiT) 
around the AFM instability in three-dimensional metals. In contrast, the power law 
(1/TT) ~ T -1 / 3 derived from the experimentally observed behaviour of (1/TT) i n a 
wide temperature range, 0.1 K < T < 10 K, was reported in [19]. 
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Figure 2. Temperature dependence of \jT{F obtained from 7 Li-NMR measurements 
on powder samples of LiV204 under different applied pressure; the data are taken 
from K. Fujiwara et al [19]. Different curves together with corresponding values of 
the fit parameter yg (0) represent the calculations based on the SCR theory. 



The discrepancy can be clearly explained by taking into account that the SCR 
theory is invoked to describe homogeneous spin fluctuations in a pure LiV204, while 
a certain amount of crystal defects and/or magnetic impurities unavoidably present in 
the measured powder samples of L1V2O4 may contribute to the nuclear spin relaxation 
as well. Actually, when approaching a magnetic instability and softening of largely 
degenerate low lying spin fluctuations, the system becomes very susceptible to weak 
perturbations including, for instance, magnetic defects. As known [13, 14, 15], the 
cooperative properties of the paramagnetic LiV 2 4 detected by NMR under ambient 
pressure can be drastically changed at sufficiently low temperatures T ~ 1 K due to a 
small amount of magnetic defects. We suggest that the low-T upturn of (1/7\T) starting 
at T upturn ~ 0.6 K is more likely to be a manifestation of the onset of spin freezing of 
magnetic defects. If so, this calls for a complete reexamination of the experimental data 
below T upturn , as done, for instance, in [15] by using a stretched exponential description 
for the nuclear spin relaxation. We expect that measurement of (l/TiT) on a sample 
with the same quality, but under the ambient pressure, would reveal the onset of spin 
freezing at a lower temperature, T upturn < 0.6 K, because of a larger distance from the 
instability. 

Finally, it is helpful to discuss the T-dependence of 1/71 predicted by the SCR 
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theory for higher temperatures, T > T . For instance, a typical behaviour 1/Ti ~ y/T is 
expected in the paramagnetic state of nearly and weakly AFM three-dimensional metals 
[21, 22]. By extending formally the present SCR theory to the region T > T , we obtain 
1/Ti ~ T a with a > 1/2, because of the SCR equation (8) reflecting an effective quasi 
one-dimensionality due to the strong anisotropy in q space of AFM spin fluctuations 
in our model. We recall, however, that the present SCR theory works well up to 40 
K, which is below T characteristic to LiV 2 4 , as evidenced from the comparison of 
theoretical results with INS data [10]. Actually, for T > T , the AFM fluctuations at 
|q| ~ Q c are suppressed and no more distinguished from those at other wave vectors in 
BZ; the system enters a spin localized regime compatible with the Curie- Weiss behavior 
of x (q = 0) observed in LiV 2 4 for T > 60 K. 

3. Summary and conclusions 

AFM spin fluctuations located in a large region of q space around the critical surface of 
a mean radius |q| ~ Q c ~ 0.6 AT 1 dominate the low-T properties of the paramagnetic 
spinel LiV 2 4 . A parametrized self-consistent renormalization theory of the AFM 
spin fluctuations was developed and applied to describe temperature and pressure 
dependences of the low-T spin-lattice relaxation rate 1/Ti i n this material. Most of 
the empirical parameters entering the present SCR theory have been estimated earlier 
and kept fixed in the present study. To simulate pressure effects in our calculations, 
the inverse static spin susceptibility x _1 {Qc) is considered to be the only fit parameter 
depending on the applied pressure. 

Comparison between NMR data and the calculated results has shown that the 
SCR theory is able to describe correctly the development of AFM spin fluctuations as 
the paramagnetic metallic state of LiV 2 4 approaches a magnetic instability under the 
applied pressure. In particular, we concluded that up to the highest applied pressure 
of 4.74 GPa, the spin system is still away from the instability. In this case, the SCR 
theory predicts the Korringa behaviour (1/TiT = constant) in a narrow temperature 
range near T = 0; on applying higher pressure, a constant value 1/TiT becomes larger 
and a temperature range near T = 0, where the Korringa relation holds, shrinks. 
The theoretical results were shown to be in a good agreement with experimental 
data for 1/TiT in wide ranges of temperature and pressure. A deviation from the 
theoretical prediction was interpreted as a signature of changing below the characteristic 
temperature T upturn (< 1 K) of the spin fluctuation dynamics. If for T > T upturn the 
nuclear spin relaxation is determined by homogeneous AFM spin fluctuations inherent 
to a magnetically pure LiV 2 4 , then below T uptU rn cooperative properties of the detected 
paramagnetic state are strongly affected by coupling of spin fluctuations to magnetic 
defects in the measured sample of LiV 2 4 . A detailed mechanism describing this 
coupling and the magnetic defect interactions in LiV 2 4 remains a challenge [14, 15] for 
a further work. 
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